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Abstract 

This library (collection of subroutines) is presented for calculating standard quantities in the de- 
composition of many-electron matrix elements in atomic structure theory. These quantities include the 
coefficients of fractional parentage, the reduced coefficients of fractional parentage as well as reduced and 
completely reduced matrix elements for several operators. So the library is assigned for any computa- 
tional scheme. The software is an implementation of a methodology based on the second quantization in 
coupled tensorial form, the angular momentum theory in 3 spaces (orbital, spin and quasispin), and the 
graphical technique of angular momentum. This implementation extends applications in atomic theory 
capabilities to partially filled /- shells and has lead to faster execution of angular integration codes. 
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PROGRAM SUMMARY 



Title of program : SQ 
Catalogue identifier: 

Program obtainable from: State Institute of Theoretical Physic and Astronomy, A. Gostauto 12, Vilnius, 2600, 
LITHUANIA. E-mail: gaigalas@itpa.lt 

Computer for which the library is designed and others on which it has been tested: 
Computers: Pentium-based PCs 

Installations: Institute of Theoretical Physics and Astronomy, A. Gostauto 12, Vilnius, 2600, LITHUANIA 

Operating systems or monitors under which the new version has been tested: Sun UNIX, LINUX 2.2.3 

Programming language used in the new version: FORTRAN 77 

Memory required to execute with typical data: 6000K Bytes 

Peripherals used: terminal, disk 

No. of bits in a word: 32 

No. of processors used: 1 

Has the code been vectorised or parallelized?: no 

No. of bytes in distributed program, including test data, etc. : 285 000 bytes 
Distribution format: compressed tar file 
CPC Program Library subprograms used: none 

Additional keywords : atomic structure, Clebsch-Gordan coefficient, configuration interaction, Reduced coefficients 
of fractional parentage, irreducible tensors, angular momentum theory in three spaces (orbit, spin and quasispin), 
second quantization in the coupled tensorial form, recoupling coefficients, complex atom, correlation, wave functions, 
bound states, LS coupling, /- shell, matrix elements, 3nj-coefficients. 

Nature of physical problem 

Accurate theoretical determination of atomic energy levels, orbitals and radiative transition data requires the calcula- 
tion of matrix elements of physical operators accounting for relativistic and correlation effects (see the multiconfigura- 
tion Hartree-Fock method [1], for example). The spin-angular integration of these matrix elements is typically based 
on standard quantities like the matrix elements of the unit tensor, the (reduced) coefficients of fractional parentage as 
well as a number of other reduced matrix elements concerning various products of electron creation and annihilation 
operators [2] . These quantities arise very frequently both in configuration interaction approaches and the derivation 
of perturbation expansions for many-particle system using symmetry-adapted configuration state functions. 

Method of solution 

This program is created involving the angular methodology of [2-7]. It has been extended to include partially filled 
/- subshells in wavefunction expansions. The classification of terms is identical to that described in [6] . 

Restrictions on the complexity of the problem 

For LS-coupling subshells states, the library provides coefficients and matrix elements for all subshells (nl) with 
Z = 0, 1, 2 and 3, and I 2 for / > 3. 

Unusual features of the program 

The library can be used as an "electronic tables" of standard quantities for evaluating general matrix elements for 
Z/S-coupled functions. 
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LONG WRITE-UP 



1 Introduction 

In order to obtain accurate values of atomic quantities it is necessary to account for relativistic and correlation 
effects (see for example [1, 2]). Relativistic effects may be taken into account as Breit-Pauli corrections or in 
a fully relativistic approach. In both cases for complex atoms and ions, a considerable part of the effort must 
be devoted to integrations over spin-angular variables, occurring in the matrix elements of the operators 
under consideration. 

Many existing codes for integrating arc based on a scheme first proposed by Fano [3]. This approach is 
based on the coefficients of fractional parentage (CFP) and then the integrations over spin-angular variables 
constitute a considerable part of the computation, especially when atoms with many open shells are treated, 
and the operators are non-trivial [4]. Over the last decade, an efficient approach for finding matrix elements 
of any one- and two-particle atomic operator between complex configurations has been developed (see papers 
Gaigalas et al [5, 6, 7, 8, 9]). It is free of the shortcomings of previous approaches (see Gaigalas [10]). This 
approach is based on a second quantization [11] and uses a coupled tensorial form for the electron creation 
and annihilation operators [12]. It also applies the theory of angular momentum [13, 14] to three different 
spaces, i.e. the space of orbital angular momentum I, spin space s and the quasispin space q [1]. The 
basic quantities of this new approach are the reduced coefficients of fractional parentage (RCFP) and the 
completely reduced matrix elements of the W^ kqklks ^ operator. 

Obviously, each computational scheme is based on a set of standard quantities to decompose the many- 
electron matrix elements. These quantities are either CFP, RCFP, the reduced matrix elements of the unit 
tensors and V^ kl \ the completely reduced matrix elements W^ kqklks, \ the completely reduced matrix 
elements of some tensorial products of second quantization operators, depending on the approach [1]. In this 
paper, we will present the library SQ (standard quantities). This library is collection of subroutines for the 
calculation of abovementioned standard quantities. The same standard quantities arise very frequently in 
effective Hamiltonian of perturbation theory or effective operator whose matrix elements between the non- 
relativistic LS coupling states are equal to the matrix elements of the full electronic hamiltonian, between the 
corresponding jj coupled relativistic states (see Rudzikas [1]), too. So the library not only supports large- 
scale computations of open-shell atoms using multiconfiguration Hartree-Fock or configuration interaction 
approaches, but may even help to develope codes for calculating the angular parts of effective operators from 
many-body perturbation theory and orthogonal operators or for evaluating relativistic hamiltonian in LS- 
coupling as well as for various versions of semi-empirical methods. The code also is intcnted for approaches 
and/or calculations presented in [15, 16, 17, 18, 19]. Some very accurate calculations was performed using 
this library (see for example [20, 21, 22, 23]), too. 

The theoretical background of this program will be presented in the section 2. It includes a brief outline 
of the quasispin concept, the defininitions of the RCFP and the reduced matrix elements of W^ kqklks \ The 
library is presented in section 3. 

2 Theoretical bacground 

The library is based on the irreducible tensorial form of the second quantization operators and on a quasispin 
technique. In this section we briefly describe this approach. 

In the quasispin representation, for a wave function of the shell of equivalent electrons \nl N aLS) a label 
Q - quasispin momentum of the shell - is introduced, which is related to the seniority quantum number 
v, namely, Q = (21 + 1 — v) /2, and its projection, Mq = (N — 21 — 1) /2. Here a denotes all additional 
quantum numbers needed for the one-to-one classification of the energy levels. Then, we can rewrite the 
wavefunction of equivalent electrons as 



\nl aQLS Mq). 
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it is possible to define the notions of a completely reduced matrix element (l aQLS\ \ \T^ k " klks ^> \ \ \l a'Q'L'S') 
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also act in 



Compared with the electron creation and annihilation operators above, the operators a 
an additional quasispin space like a tensor component with rank q and a projection m q = There is the 
following relation known between the reduced matrix element of a creation operator and the CFP [24]: 



(l N aQLSWa^W^- 1 a'Q'L'S 1 ) = {-l) N y/N [J,L,S\{l N aQLS\\l N ~ l (a'Q'L'S') l) 
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where [L, S] = (2L + 1) (2S + 1). Eqs. (4) and (2) can be used to define the relation between the CFP and 
its reduced counterpart in Q-space. Introducing the z— projection, Mq, of the quasispin, this relation is 
given by [1] 
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Tables of numerical values of (l aQLS\\\a^ qls ^\\\l aQ'L'S') are presented in Rudzikas and Kaniauskas [12] 
when I = 0,1,2. For the tensorial product of two one-electron operators, the submatrix element equals 
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On the right-hand side of equations (5) and (6) only the Clebsch-Gordan coefficient 
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{nl aQLS\\\W^ klk ^\\\nl a'Q'L'S') denotes reduced in quasispin space submatrix element (completely 

reduced matrix element) of the triple tensor W^ klk2 ^ — [a^ 1 ^ x a^ qls ^ 
a following way: 
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So, by applying the quasispin method for calculating the matrix elements of any operator, we can use 
the RCFP or the completely reduced matrix elements of W^ klk ^ (nl,nl), which arc independent of the 
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occupation number of the shell. The main advantage of this approach is that the standard data tables in 
such a case will be much smaller in comparison with tables of the usual coefficients and, therefore, many 
summations will be less time-consuming. Also one can see that in such an approach the completely reduced 
matrix elements of standard tensors and RCFP actually can be treated in a uniform way as they all arc 
the completely reduced matrix elements of the second quantization operators. Hence, all methodology of 
calculation of matrix elements will be much more universal in comparison with the traditional one. 



3 Description of the Library 
3.1 SAI_SQLS1 

The section SQLS1, (standard quantities in LS coupling, part one) is a collection of utilities for calculation 
of standard quantities as: 
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• The Qj— and 9j— coefficients. 

• The RCFP ( I QLS ||| a^ s ) ||| I Q'L'S 1 ). 

• Reduced matrix element ( I QLS \\\ [ a (qls ^ x ]( k ^ 3 ) ||| t Q'L'S' ). 

• The matrix elements of type ( l N QLS \\ [a%$ x a& l £]V* k -)\\ l N ' Q'L'S' ), 
( l N QLS || [ati$ x [ x || ^ Q' L ' S > ), 

( l N QLS || [[ x ]^ x ](*<*.) \\ l N ' Q'L'S' ), 
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The subroutines presented in this section are independent and may be useful for other programs. 67 
subroutines are contained in this library. Was created the library [25] on the low-powered computers BESM- 
6 of the 1980s for some similar quantities for the s-, p- and d- shells. It was written in Fortran 4 and had 
17 subroutines. It was based on the approach [26] and assign for calculation submatrix elements of the 
nonrelativistic effective Hamiultonian of atom in the first two orders of stationary perturbation theory in 
the case of configuration of two open-shells nil^ 1 ri2^ 2 . The present library is more general and has more 
future because it is useful for different approaches with any number of open shells with 1 = 0, 1, 2 and 3, 
and I 2 for I > 3. 



3.1.1 Common Blocks 

Most of the subroutines use common block GLCONS from MCHF atomic structure package [2, 27]. The 
parameters contained in this block are defined in BLOCK DATA GLCONS. 

The Common blocks /MT/ and /MT67/ are important, too. The first one has the array MT(40), which 
contains all the term characteristics of s-, p- and d- shells, that are needed while calculating matrix elements 
in the quasispin formalism (see Table 1). An element of the MT array is indicated in the column No of Table 
1, and the Term column indicates the characteristics of a term contained in that element of MT array. In 
other words, for all the s-, p- and d- shells, the terms are numbered from 1 to 40. The terms are marked as 

(2S+1) £ (20+1) in the taWe ^ 

All the terms for /- shell (see Table 2) are similarly placed in the Common block /MT67/. The No 
column of Table 2 indicates the term number. The TERM column has ( 2S + 1 )i J Nr i n it, the 2Q column has 
the quasispin momentum Q multiplied by two. See [8, 28] for the details of the classification of /- shell 
terms. 

The Common block /MT/ is defined in BLOCK DATA TERMLS. Meanwhile the /MT67/ is defined 
in BLOCK DATA TRMF from the library SAI_SQLS2 (see below). BLOCK DATA TERMLS defines 
COMMON /SKMT2/, which contains term characteristics for special cases, i.e. for shells I — 3 — 9 and shell 
occupation numbers N = 1,2. 

Single shell data are stored in the two arrays I and B. The former consists of 
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Table 1: Allowed couplings of l N states for I = — 2. The subshell quasispin angular momentum Q, spin 
angular momentum S, the subshell angular momentum L and seniority quantum number v are denoted 
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• 1(1) is the state number of the shell (see Tables 1,2). 

• 1(2) is the principal quantum number n. 

• 1(3) is the orbital quantum number I. 

• 1(4) is the number of electrons in the subshell. 

• 1(5) is the shell total angular momentum L multiplied by two. 

• 1(6) is the shell total angular momentum S multiplied by two. 

• 1(7) is the shell total quasispin Q multiplied by two. 
The array B contains 

• B(l) is the shell quasispin Q. 

• B(2) is the shell total angular momentum S. 

• B(3) is the shell quasispin projection Mq. 

These arrays are placed in Common blocks /TRK/ and /TRK2/. In particular, these are: 
Name Dimension Function 



/TRK/ The data of the orbitals for first two shells 

BDS1 3 the array B for the first shell of the ket function 

BDS2 3 the array B for the second shell of the ket function 

BKS1 3 the array B for the first shell of the bra function 

BKS2 3 the array B for the second shell of the bra function 

IBDS1 7 the array I for the first shell of the ket function 

IBDS2 7 the array I for the second shell of the ket function 

IBKS1 7 the array I for the first shell of the bra function 

IBKS2 7 the array I for the second shell of the bra function 

/TRK2/ The data of the orbitals for last two shells 

BDS3 3 the array B for the third shell of the ket function 

BDS4 3 the array B for the fourth shell of the ket function 

BKS3 3 the array B for the third shell of the bra function 

BKS4 3 the array B for the fourth shell of the bra function 

IBDS3 7 the array I for the third shell of the ket function 
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Table 2: Allowed couplings of f N . The subshell spin angular momentum S, number Nr and the subshcll 
angular momentum L are denoted ( 2S + 1 )£ Ar? \ The quasispin angular momentum is denoted Q. 
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Table 2 (continued) 
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IBDS4 7 the array I for the fourth shell of the kct function 

IBKS3 7 the array I for the third shell of the bra function 

IBKS4 7 the array I for the fourth shell of the bra function 



In this library the common blocks /MEDEFN/ and /FACT/ from the earlier version of the MCHF 
atomic-structure package [2, 27], and the newly created auxiliary COMMON blocks /KAMPAS/, /PJBOLS/, 
/RIBOLSF/, /PJBOF/, /PJBOLS3/, are used. 

The authors would like to draw attention at the ordering of terms in Table 2, which is tuned to simplifying 
the placement of tables of /- shell reduced coefficients of fractional parentage into DATA blocks. 



3.1.2 Subroutines 

The subroutine NUMTER 

FUNCTION NUMTER has several modes of operation, but only one of these is possible while using it in 
other programs. If the input values are the shell's total quasispin Q multiplied by two (input argument 
I2Q), the shell's total spin S multiplied by two (input argument I2S), the shell's total angular momentum L 
multiplied by two (input argument I2L) and the quantum number I, it will find the number of the s-, p- d- 
shell's term, as numbered in Table 1. The other arguments of this subroutine, NK and ND, should be set 
to 3. In this mode, a COMMON block /MT/, is needed by the program, which is defined in BLOCK DATA 
TERMLS. 



The subroutine RUMT 

The subroutine RUMT has several modes of operation, but only one of these is possible while using it in other 
programs. For an input of the orbital quantum number I (input argument LL) and the term number from 
Table 1 (input argument KNT), it finds the shell's total quasispin Q multiplied by two (output argument 
LQ), the shell's total spin S multiplied by two (output argument LS) and the shell's total angular momentum 
L multiplied by two (output argument L). It finds these characteristics for s—, p- and d— shells. A COMMON 
blocks /MT/ is needed by this program, which is defined in BLOCK DATA TERMLS. 

The subroutine C0T5S 

This subroutine determines the value of the Clcbsch Gordan coefficients: 



Q 



c 



QM SM CM 



(8) 



The subroutine has the input arguments Q, QM, SM, C, CM and output argument A. The subroutine 
performs its calculations by employing analytical expressions from Varshalovich et al [29] . 

The subroutine C1E0SM 

This routine determines the value of the Clebsch-Gordan coefficients: 



Q 1 C 
QM CM 



(9) 



The subroutine has the input arguments Q, QM, C, CM and output argument A. The subroutine performs 
its calculations by employing analytical expressions from Varshalovich et al [29] . 



The subroutine C1E1SM 

This subroutine determines the value of the Clebsch-Gordan coefficients: 



Q 

QM 



C 
CM 



(10) 



The subroutine has the input arguments Q, QM, C, CM and output argument A. The subroutine performs 
its calculations by employing analytical expressions from Varshalovich et al [29] . 
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The subroutine SIXJ 



This routine determines the value of the 6 j -coefficients: 

f 1/2 J/2 K/2 \ 

\ L/2 At/2 N/2 J' [ ! 

The subroutine has the input arguments /, J, K, L, M , N , ITIK and output argument 5*/. If the parameter 
ITIK = 0, the subroutine does not check the triangular conditions for 6 j -coefficient. In other cases it checks 
these. If any of the parameters of 6j-coefficient is equal to 0, \, 1, §, 2, 3, 4, the subroutine calculates the 6j- 
cocfficients according to analytical formulas [29, 14]. Otherwise, the customary calculations are performed. 
In that case, the COMMON block /FACT/ must be defined. This is done by addressing the SUBROUTINE 
FACTRL from the library MCHF_LIB_COM [2, 27]. 



The subroutine NINE 

Jl/2 J2/2 J3/2 ] 

Ll/2 L2/2 L3/2 \. (12) 
Kl/2 K2/2 K3/2 J 

The subroutine has the input arguments Jl, J2, J3, LI, L2, L3, Kl, K2, K3, I and output arguments IN 
and A A. If the parameter 1=1, the subroutine only checks the triangular conditions of a 9j- coefficient. If 
these are not satisfied, then IN = 0, and IN = 1 otherwise. At other values of /, the subroutine calculates 
the value of 9j- coefficient and assigns it to the output parameter AA. 

The subroutine SLS 

This subroutine determines the value of the RCFP: 

( I QLS HI a (9/s) HI I Q'L'S' ) (13) 

for p-, d-, and /- shells (see (39) in [6]). The routine uses the table of reduced matrix elements of the a^ qls ^ 
tensor operator from [8]). The subroutine has the following arguments: 

1. L is the orbital quantum number I. 

2. IT is the state number of the bra function (see Tables 1,2). 

3. LQ is the quasispin Q for the bra function multiplied by two. 

4. LL is the total angular momentum L for the bra function multiplied by two. 

5. LS is the total angular momentum S for the bra function multiplied by two. 

6. ITS is the state number of the ket function. 

7. LQS is the quasispin Q for the ket function multiplied by two. 

8. LLS is the total angular momentum L for the ket function multiplied by two. 

9. LSS is the total angular momentum S for the ket function multiplied by two. 

10. S is the value of the reduced matrix element (13) which is returned by the subroutine. 



The subroutine RWLS 

The routine determines the value of the reduced matrix element: 

( I QLS HI [ a (9 ' s) x a {qls ^ }( k ^ k ^ \\\ I Q'L'S' ). (14) 

The routine uses the tables of reduced matrix elements of the tensor operator for s-, 

p- and d— subshells (see Spakauskas et al [30]), and for the /- shell the expression (34) from paper [5] is 
used. The subroutine does not calculate the simple case of k\ — fc 2 = k 3 = 0, because then the operator is 
just [ x ](° °) = -{21 + l) 1 / 2 (expression (15.54) in Rudzikas [1]). The subroutine has the formal 

arguments: 
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1. Kl is the rank k\. 

2. K2 is the rank fc 2 . 

3. K3 is the rank k 3 . 

4. L is the orbital quantum number I. 

5. Jl is the state number of the bra function (see Tables 1,2). 

6. J2 is the state number of the ket function. 

7. W is the value of the reduced matrix element (14) which is returned by the subroutine. 
The subroutine Wl 

This subroutine determines the value of the matrix element: 

( l N QLS || [<f x <'>]( fc ' fc ->|| l N ' Q'L'S 1 ). (15) 

While calculating cases where the orbital number 1=0, 1, 2, 3 and the shell occupation number N > 2, the 
program relics on the expression (31) from the paper [5]. In that case, the subroutine finds the Clebsch- 
Gordan coefficient which gives the dependence on the shell occupation number. If the tensor product (15) 
consists of either two electron creation or two annihilation operators then C1E1SM is called. Otherwise 
CLEOSM is called. The subroutine RWLS finds the reduced matrix elements of the operator 
a (qls) ](fe;*=s)^ j n other cases, the program calculates according to the expression (40) from paper [6]. The 
subroutine has the formal arguments: 

1. IK is the array I for the bra function. 

2. BK is the array B for the bra function. 

3. ID is the array I for the ket function. 

4. BD is the array B for the ket function. 

5. K2 is the rank h\. 

6. K3 is the rank k s . 

7. QM1, QM2 are the quasispin projections in (15). 

8. W is the value of the reduced matrix element (15) which is returned by the subroutine. 
The subroutine AWP1LS 

The routine determines the value of the matrix elements: 

( l N QLS || [a&W x [ a<£? x ](*!*)](*«*.) || l N ' Q'L'S' ). (16) 

While calculating cases where the orbital number 1=0, 1, 2, 3 and the shell's occupation number N > 2, 
the program relies on the expression (31) from the paper [5]. In that case, the subroutine IZAS1 checks 
that the subshell has a state with the specified characteristics. The subroutine ITLS2 finds the first and 
the last numbers of the state from the running intermediate sum in array MT. RUMT finds the shell's 
total angular momentum LS and quasispin Q for each intermediate state. The routine IXJTIK checks all 
triads. The subroutine C0T5S finds the Clebsch-Gordan coefficient giving the dependence on the shell's 
occupation number and SLS finds the reduced matrix element of a( qls *> tensor operator (see in Section 3.1). 
The second part of the expression is calculated by the routine Wl (see in Section 3.1.2). The routine SIXJ 
finds 6j-symbol. In other cases, the program calculates according to the expression (40) from paper [6]. The 
subroutine has the arguments: 

1. IK is the array I for the bra function. 

2. BK is the array B for the bra function. 
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3. ID is the array I for the ket function. 

4. BD is the array B for the ket function. 

5. Kl is the rank k\. 

6. K2 is the rank fc 2 . 

7. K3 is the rank ki. 

8. BK4 is the rank k s . 

9. QM1, QM2 and QM3 are the quasispin projections in (16). 

10. AW is the value of the reduced matrix element (16) which is returned by the subroutine. 
The subroutine WAP1LS 

This subroutine determines the value of the matrix elements: 

( l N QLS || [[ <f x P 1 ^ x }^ || l N ' Q'L'S' ). (17) 

The structure of the routine WAP1LS is the same as that of AWP1LS. The subroutine has the formal 
arguments: 

1. IK is the array I for the bra function. 

2. BK is the array B for the bra function. 

3. ID is the array I for the ket function. 

4. BD is the array B for the ket function. 

5. Kl is the rank k\. 

6. K2 is the rank fc 2 . 

7. K3 is the rank fc;. 

8. BK4 is the rank k s . 

9. QM1, QM2 and QM3 are the quasispin projections in (17). 

10. WA is the value of the reduced matrix element (17) which is returned by the subroutine. 
The subroutine WWLS1 

The routine determines the value of matrix elements: 

( l N QLS || [ [ <f x a*) ](*'*•> x [ x ]( fe ^)](°°) || l N ' QLS ). (18) 

The routine WWLS1 uses the routines ITLS, RUMT, IZAS1 for calculation of this sort of the matrix element. 
The subroutine Wl calculates the first and the second parts of the operator calculates. The subroutine has 
the formal arguments: 

1. IK is the array I for the bra function. 

2. BK is the array B for the bra function. 

3. ID is the array I for the ket function. 

4. BD is the array B for the ket function. 

5. K2 is the rank fc;. 

6. K3 is the rank k s . 

7. QM1, QM2, QM3 and QM4 are the quasispin projections in (18). 

8. WW is the value of the reduced matrix element (18) which is returned by the subroutine. 
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The subroutine WWPLS1 



The routine determines the value of matrix elements: 

( l N QLS || [ [ x ]< fclfa > x [ x ](^)](Ws) || Z"' QLS ). (19) 

The routine WWPLS1 uses the routines ITLS, RUMT, IZAS1 for calculation of this sort of the matrix 
element. The subroutine Wl calculates the first and the second parts of the operator calculates. The routine 
SIXJ finds 6j-symbol. The subroutine uses the COMMON block /TRK/, which contains arrays I and B for 
the bra and ket function (ID1, IK1,BD1 and BK1). It uses the expression (40) from [6]. The subroutine has 
the formal arguments: 

1. Kl is the rank k\. 

2. K2 is the rank ki. 

3. K3 is the rank /c 3 . 

4. K4 is the rank k^. 

5. K5 is the rank fcj. 

6. K6 is the rank k s . 

7. QM1, QM2, QM3 and QM4 are the quasispin projections in (19). 

8. WW is the value of the reduced matrix element (19) which is returned by the subroutine. 

3.2 SAI_SQLS2 

The library SAI_SQLS2 is the second part of Standard quantities in LS- coupling. Everything in it is related 
to the /- shells. The tables of CFP for /- shell (see [8]) and the term characteristics are in this library. This 
library can be used in other programs to its full extent, by employing programs from SAI_SQLS1 library. 
These are the subroutines SLS, RWLS, Wl, AWP1LS, WAP1LS, WWLS1, WWPLS1. The remaining two 
subroutines may also be used independently. 

The subroutine NUMTERF 

The subroutine has several modes of operation. But only one of these is useful while operated independently 
from the MCHF atomic structure package [2, 27]. In that case, after giving number Nr from Table 2 (input 
argument I2N), the values of the shell total spin S multiplied by two (input arguments I2S and N), the 
shell total angular momentum L multiplied by two, (input argument I2L) and the shell total quasispin Q 
multiplied by two (input argument I2Q) for the FUNCTION NUMBER, it finds the number of /- shell term, 
as numbered in Table 2. The subroutine needs the COMMON block /MT67/, which is defined by BLOCK 
DATA TERMF. 

The subroutine RUMT67 

For an input of the term member from Table 2 (input argument KNT), the subroutine RUMT67 finds the 
number Nr from Table 2 (output argument NR), the shell total quasispin Q multiplied by two (output 
argument LQ), the shell total spin S multiplied by two (output argument LS) and the shell total angular 
momentum L multiplied by two (output argument L). The subroutine needs the COMMON block /MT67/ 
which is defined by BLOCK DATA TERMF. 

4 Examples 

A driver program (see Figure 1) accompanying the library SQLSF illustrates three examples. 

The examples show that some of the subroutines contained in the libraries may serve as an electronic 
version of Nielson and Koster [31] or [8] tables. As such they may serve as a basis for extending the capabilities 
of programs that rely on the principle of calculating the determinants to arbitrarily filled /- shells (see for 
example Eissner et al [32] or Zatsarinny [33]). 

Here we describe each example separately, the EXAMPLES RUN OUTPUT being shown in Table 4. 
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IFCICASE .Eg. 1) THEN 

WRITEC6 , ' (/A) ' ) ' orbital quantum number 1 (II) ' 
READ (5, ' (ID ') L 
IF(L.LT.3) THEN 

WRITEC6 , ' (/A) ' ) ' 2*Q 2*L 2*S for bra function (312) ' 

READ (5, ' (312) ') IQB.ILB.ISB 

JB=NUMTER ( I QB , I SB , I LB , L , 3 , 3) 

WRITEC6 , ' (/A) ' ) ' 2*Q 2*L 2*S for ket function (312) " 
READ (5, ' (312) ') IQK, ILK, ISK 
JK=NUMTER(iqK,ISK,ILK,L,3,3) 
ELSEIF(L . EQ . 3) THEN 

WRITEC6 , ' (/A) ' ) ' 2*Q 2*L 2*S Nr for bra function (412) ' 
READ (5, ' (412) ') IQB , ILB , ISB , INB 
JB=NUMTERF (INB , ISB , ILB , ISB , IQB) 

WRITEC6 , ' (/A) ' ) ' 2*q 2*L 2*S Nr for ket function (412) ' 
READ (5, ' (412) ') iqK, ILK, ISK, INK 
JK=NUMTERF (INK , ISK , ILK , ISK , IQK) 
ENDIF 

CALL SLS (L , JB , IQB , ILB , ISB , JK , IQK , ILK , ISK , S) 
WRITE (6, ' (/A,F17.7) ') ' Value= \S 
ELSEIFdCASE .EQ. 2) THEN 

WRITEC6 , ' (/A) ' ) ' orbital quantum number 1 (II) ' 
READ (5, ' (II) ') L 

WRITEC6 , ' (/A) ' ) ' number electrons in the shell N (II) ' 
READ (5, ' (II) ') N 

WRITEC6 , ' (/A) ' ) ' bra and ket number of term (213) ' 
READ (5, ' (213) ') JB,JK 
IF(L.LT.3) THEN 

CALL RUMT(JB,L,IQB,ISB,ILB) 

CALL RUMT(JK,L,iqK,ISK,ILK) 
ELSEIF(L . EQ . 3) THEN 

CALL RUMT67(JB,NRB,iqB,ISB,ILB) 

CALL RUMT67(JK,NRK,IQK,ISK,ILK) 
ENDIF 

CALL SLS (L , JB , IQB , ILB , ISB , JK , iqK , ILK , ISK , S) 

QK=DBLE(iqK)*0.5 

QB=DBLE(iqB)*0.5 

QMK=DBLE(N-1-2*L-1) *0 . 5 

QMB=DBLE(N-2*L-1) *0 . 5 

QM=0 . 5 

CALL C0T5S (QK , QMK , qM , QB , QMB , A) 

CFP=S*A/DSQRT(DBLE(N* (IQB+1) * (ILB+1) * (ISB+1) ) ) 
IF(((N-1)/2)*2.NE.N-1) CFP=-CFP 
WRITE(6, ' (/A.F17.7) ') ' Value= > ,CFP 
ELSEIFdCASE .EQ. 3) THEN 

WRITEC6 , ' (/A) ' ) ' orbital quantum number 1 (II) ' 
READ (5, ' (ID ') L 

WRITEC6 , ' (/A) ' ) ' number electrons in the shell N (ID ' 
READ (5, ' (ID ') N 

WRITEC6 , ' (/A) ' ) ' first and last number of term for bra (213) 
READ (5, ' (213) ') JFB , JLB 

WRITEC6 , ' (/A) ' ) ' first and last number of term for ket (213) 
READ (5, ' (213) ') JFK, JLK 

WRITE(6 , ' (/A) ' ) ' the ranks Kl K2 K3 (311) ' 
READ(5, ' (311) ') K1,K2,K3 
WRITE(6, ' (311) ') K1,K2,K3 
DO 2 JTB=JFB , JLB 
IF(L.LT.3) THEN 

CALL RUMT(JTB,L,IQB,ISB,ILB) 
ELSEIF(L . EQ . 3) THEN 

CALL RUMT67(JTB,NRB,IQB,ISB > ILB) 
ENDIF 

DO 3 JTK=JFK, JLK 
IF(L.LT.3) THEN 

CALL RUMT(JTK,L,IQK,ISK,ILK) 
ELSEIF(L.EQ.3) THEN 

CALL RUMT67(JTK,NRK,IQK,ISK,ILK) 
ENDIF 

CALL RWLS(K1,K2,K3,L,JTB,JTK,W) 

QK=DBLE(IQK)*0.5 

QB=DBLE(IQB)*0.5 

QM=DBLE(N-2*L-1) *0 . 5 

CALL C1E0SM(QK,QM,QB,QM,A) 

VK=-0 . 5*W*A/DSQRT(DBLE( (IQB+1) * (2*K2+1) ) ) 

WRITE(6, ' (/A,/2I3,F17.7) ') 
: ' number for bra : number for ket : value ',JTB,JTK,VK 

3 CONTINUE 
2 CONTINUE 
ENDIF 



Figure 1: The program for three examples 
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Case 1. Finding the RCFP 

In this case the subroutine finds the value of a RCFP ( / 8 5° ||| ||| / 7 F° ). From Table IV of [8] we 

see that the value of this coefficient is equal to 

( / 8 5° HI HI / 7 F° ) - -4V7. 

This is in accordance with the result obtained with SLS subroutine (see in Section 3.1). 

It is necessary to mention that in the calls to this and some other subroutine, the input of a term number 
and term characteristics is needed. In this case the user indicates only the term characteristics, and finds 
the term number by using the subroutine NUMTER (see in Section 3.1) or NUMTERF (see in Section 3.2). 
In other examples we will present a way, showing how to call subroutines of analogous type, by defining just 
the term number by the user. 

Case 2. Calculating the CFP 

This case illustrates the finding of a CFP value, using the SLS program. It may be useful for programs that 
are based on the tables of CFP, which are much more extended than the tables of RCFP. The values of CFP 
are found by programming the expression (5). 

The value of a CFP ( f 6 D° || f 6 5 P° f ) is found in the example. We see from Nielson and Kostcr [31] 
that the value of this CFP is equal to 

{f 6 ^°||/ 6 5 P° f) = ~. 

This is in accordance with the result obtained with SLS subroutine. 

See [8, 28] for details of the encoding of the /- shell terms. See [8] for a more extensive use of the tables 
of RCFP in finding the CFP. In that paper, also the problems related to that task are discussed and the 
ways to solve these are indicated. 

Case 3. Calculating the matrix elements of V k operator 

In this section it is demonstrated how to calculate the V 11 operator, using the subroutine RWLS (see in 
Section 3.1), which calculates the reduced matrix elements (14). For that purpose the relation (16.34) from 
Rudzikas [1] is used. While using the U k and V kl tables, one must pay attention to various phase conventions 
used in the literature. In addition, small differences in the definitions of U k occur. Some authors, Karazija 
et al [34] among them, tabulate the submatrix elements 

(l N aSL\\U k \\l N a'S'L'), (20) 

others, like Nielson and Koster [31] or Cowan [35], tabulate 

{l N aL\\U k \\l N a'L'), ( 21 ) 

although they use the notation of (20). Meanwhile, the relation between these two coefficients is: 

(l N aSL\ \U k \ \l N a'S'L') = S(S, S')y/(2S + l)(l N aL\\U k | \l N a'L'). (22) 

The submatrix elements are defined as (20) if we use the relations between matrix elements of W^ klk2k3 ^ and 
U k as presented in the Rudzikas monograph [1]. Reduced matrix elements of the operator V 11 (/ 7 6 P°\ l^ 11 1 1/ 7 4 S 1 ) 
are calculated in the example. The numerical value of this reduced matrix element is taken from the Nielson 
and Koster [31] tables 

(/ 76 P°||T/ 11 || / 74 ^) =-yjl 

and agree with our value. 
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Table 4. EXAMPLES RUN OUTPUT 

Case 1. Finding a reduced coefficient of fractional parentage 

» Example 



EXAMPLES for LIBRARIES 



Which case? (1,2,3) 
>1 
CASE 1 



orbital quantum number 1 (ID 
>3 

2*Q 2*L 2*S Nr for bra function (412) 

> 7 

2*Q 2*L 2*S Nr for ket function (412) 

> 1 6 6 

bra (number 1 2*Q 2*L 2*S) 

1 3 7 
ket (number 1 2*Q 2*L 2*S) 

120 3 1 6 6 
Value= -10.583005244258363 
END OF CASE 1 



Case 2. Calculation of the coefficient of fractional parentage 

>> Example 



EXAMPLES for LIBRARIES 



Which case? (1,2,3) 
>2 
CASE 2 



orbital quantum number 1 (ID 
>3 

number electrons in the shell N (ID 
>7 

bra and ket number of term (213) 
3129 

bra (number 1 2*Q 2*L 2*S N) 
3 3 4 5 7 
ket (number 1 2*0. 2*L 2*S) 

129 3 1 2 4 
Value= -0.247435829652697 
END OF CASE 2 



Case 3. Calculation of the matrix elements of operator V k 

» Example 



EXAMPLES for LIBRARIES 



Which case? (1,2,3) 
>3 
CASE 3 



orbital quantum number 1 (ID 
>3 

number electrons in the shell N (II) 
>7 

first and last number of term for bra (213) 

> 2 2 

first and last number of term for ket (213) 

> 8 8 

the ranks Kl K2 K3 (311) 
>111 
111 

number for bra : number for ket : value 

2 8 -0.534522483824849 
END OF CASE 3 
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5 Conclusion 



Accurate theoretical determination of atomic energy levels, orbitals and radiative transition data requires the 
calculation of matrix elements of physical operators accounting for relativistic and correlation effects. The 
spin-angular integration of these matrix elements is typically based on standard quantities like the matrix 
elements of the unit tensor, the (reduced) coefficients of fractional parentage as well as a number of other 
reduced matrix elements concerning various products of electron creation and annihilation operators. These 
quantities arise very frequently in large-scale computations of open-shell atoms using multiconfiguration 
Hartree-Fock or configuration interaction approaches, in calculating the angular parts of effective operators 
from many-body perturbation theory or for evaluating relativistic hamiltonian in LS-coupling as well as 
in various versions of scmiempyrical methods. The library SQ is assigned for the calculation of all these 
standard quantities. It can be used as "electronic tables" of standard quantities for evaluating general 
matrix elements for LS'-coupled functions, too. 

The library SQ is created involving the angular methodology of [5, 6, 7, 8, 9, 10]. For LS'-coupling 
subshells states, the library provides coefficients and matrix elements for all subshells (nl) with I = 0, 1, 2 
and 3, and I 2 for I > 3. 

Program is obtainable from State Institute of Theoretical Physic and Astronomy, A. Gostauto 12, Vilnius, 
2600, LITHUANIA. E-mail: gaigalas@itpa.lt. 
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"Prongraminc biblioteka atomy strukturos teorijos standartiniams dydziams skaiciuoti." 
G. Gaigalas 

Santrauka 

Straipsnyje pateikta biblioteka (paprogramiy rinkinys), parasyta Fortran 77 programavimo kalba. Ji 
skirta standartiniy dydziy skaiciavimui, su kuriais susiduriame atomo teorijoje nagrinejant fizikiniy bei 
efcktyviniu operatoriy matricinius elementus LS rysyje. Visu pirma, tai kilminiai koeficientai bei viene- 
tiniy tenzoriy U k ir V kl matriciniai elementai, kurie paprastai pasirodo nagrinejant matricinius elementus 
klasikiniais metodais. Be to biblioteka skirta ir tokiu standartiniy dydziy nagrincjimui, kurie pasirodo 
naudojant kvazisukinio formalizma. Tai subkilminiai keoficientai, bei redukuotiniai kvazisukinio erdveje, 
vicnctiniai tcnzoriai W^ kqk ' ks \ Taigi darbe pateikta biblioteka yra bendro pobudzio ir skirta matriciniy 
elementy skaiciavimui bet kokia mokslineje literaturoje zinoma mctodika. Pati biblioteka yra parasyta 
remiantis antrinio kvantavimo suristamc tcnzoriniamc pavidalc formalizmu, judcsio kickio momcnto tcorija 
trijose erdvese (orbitiniamc, sukininiame ir kvazisukininiame) ir grafine judesio kiekio momento teorija. Si 
biblioteka isplecia atomo teorijos galimybes, kadangi joje yra pilnai itraukti f sluoksniai su bet kokiu sluoksnio 
uzpildymo skaiciumi, bei leidzia pagreitinti jau esamas programas. 
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